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We consider the problem of deciding if a given three-party entangled pure state can be converted, 
with a non-zero success probability, into a given two-party pure state through local quantum op- 
erations and classical communication. We show that this question is equivalent to the well-known 
computational problem of deciding if a multivariate polynomial is identically zero. Efflcient ran- 
domized algorithms developed to study the latter can thus be applied to the question of tripartite 
to bipartite entanglement transformations. 
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A basic question concerning quantum entanglement is 
whether it can be transformed in a particular manner. 
More specifically, given a multipartite system originally 
in some entangled state, what are the possible final states 
the system can realize if the only allowed quantum op- 
erations are performed locally, one subsystem at a time 
while assisted by global classical communication? Such 
protocols are called Local Operations and Classical Com- 
munication (LOCC), and any LOCC transformation that 
yields a certain final state with just a nonzero probabil- 
ity is called Stochastic LOCC (SLOCC). For some initial 
state IV') and target state |(/)), an SLOCC conversion be- 
tween \ij}) and 10) is denoted as |i/') ^^^^^ |0). In this 
Letter, we study the computational problem of deciding 
if |-)/i) 10) ig feasible, given the classical description 

of It/i) and 10). 

In general, deciding the feasibility of state transfor- 
mations is difficult, with the challenge becoming more 
formidable as more parties are involved. This increase 
in difficulty can be made precise by using the language 
of computational complexity theory which groups 
problems according to the amount of resources needed 
to solve them. For SLOCC transformations of bipartite 
pure states, the computational resources in deciding con- 
version feasibility increases polynomially in the dimen- 
sion of either subsystem undergoing the transformation 
thus the SLOCC bipartite conversion problem be- 
longs to the complexity class P. On the other hand, it 
was recently shown that the SLOCC convertibility be- 
tween tripartite states is an NP-Hard problem, implying 
that no polynomial time decision algorithm exists if 
NP, where NP is a central complexity class that includes 
numerous naturally occurring problems not known to be 
in P. These results can be interpreted as a formal indica- 
tion that it is generally much easier to study pure state 



entanglement transformations in bipartite systems than 
in tripartite systems. 

A natural next step is to examine a hybridization 
of these two transformation classes and see where the 
problem of tripartite to bipartite SLOCC entanglement 
conversions fits in the complexity spectrum. This is a 
subset of tripartite transformations in which Alice, Bob 
and Charlie initially share a three-way entangled state 
\^)abc^ but they end with a state \4')ab where Alice 
and Bob arc still entangled but Charlie is completely un- 
entangled from the other two. 

In this Letter, we show that deciding tripartite to bi- 
partite convertibility is equivalent to Polynomial Identity 
Testing (PIT), which is the task of determining whether 
two polynomials given in algebraic formulas are equiva- 
lent (or cquivalently, if a polynomial is identically zero) . 
PIT is a classical problem in theoretical computer sci- 
ence with many important applications, such as in per- 
fect matching 01, multiset equality testing 0], and pri- 
mality testing [5|. In particular, it is known that PIT 
admits a polynomial time randomized algorithm but is 
not known to have a deterministic polynomial time al- 
gorithm. The failure to "derandomize" the algorithm is 
shown to arise from the difficulty of proving super- 
polynomial lower bounds on general computation mod- 
els: if PIT is in P. then some other problems would not 
have an efficient algorithm. Thus PIT has been a central 
problem to the fundamental (and open) question if ran- 
domness is useful in computation (i.e. the BPP versus 
P question) . The equivalence of the convertibility ques- 
tion to PIT implies the former also admits a randomized 
polynomial time algorithm, while it will remain a diffi- 
cult question if it has a polynomial time algorithm. To- 
gether with the previously mentioned complexity results 
on other convertibility questions, those findings establish 
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close relationships between quantum entanglement and 
computational complexity theory. 

In the rest of this Letter, we present the formal state- 
ment of our result and its proof. We start with the nec- 
essary notation and some observations. For a tripartite 
pure state \'4')abCi Alice and Bob's joint state is de- 
scribed by = Trc{ABc\tp){'tp\ABc), with a mixed 
state representation X]i=i-Pi|6i)(ei|i where {ei\ej) = 5ij 
and Pi > 0. The "subnormalized" eigenstates {\ei) = 
^ypi\ei)}i=i...n span the space supp{p^g) called the sup- 
port of Pj^g- An ensemble of pure states {\qj)}j=i...t 
satisfies X)i=i ~ Pab ^^'^ only if there exists a 



unitary matrix U 



Hj\l<i,j<t, 



such that 



1=1 



(1) 



where \ei) is the zero vector for i > n Furthermore, 
there is a one-to-one correspondence between a measure- 
ment of Charlie consisting of rank-one measure operators 
and a unitary matrix U that defines the resultant pure 
state ensemble through Eq. ([T|) shared by Alice and Bob 
following Charlie's measurement. 

A consequence of these facts is contained in the fol- 
lowing theorem whose statement requires one final bit of 
terminology. For a bipartite pure state |(^), its Schmidt 
rank refers to the minimum number of product states 
needed to express it and is denoted by rk{\(j))). Equiva- 
lently, for arbitrary bases {\i)A}i=i---dA and {|i)s}i=i - dB 
of Ha and Hb respectively, \(p) can be uniquely identi- 
fied with a cIb ^ d-A matrix 4> mapping Ha to Hb by 
10) = (1(8)$) \i)A\i)A, and rk{\(f>)) equals the ma- 

trix rank of <f>. 



Theorem 1. Let 

Schmidt rank d. 



(j)) AB be a bipartite pure state with 
Then for tripartite state \'4') ABC , 
VP) ABC \<t')AB if 0-Tid only if there exists \(j)')AB G 

supp(pyis) such that rk{\(j}') ab) > d- 



Proof. Let \iI})abc = J27=i\^i)AB\ei}c be a Schmidt 
decomposition of \iIj)abc- (=^) If the transformation is 
possible, then there are linear operators A, B, C such 
that A(g) B (g) C\iP)abc = \<I>)ab\0)c 01 • Consequently, 
\cj,)AB = A® for 10') = EtiiCh) c)\e^)AB- 

Since \(f>') e supp(p^B) and A iSi B cannot increase the 
Schmidt rank of \4>'), we have rk{\(f)')) > d. (<^) Con- 
versely, assume the existence of |0') £ supp(pAs) with 
(0'|0') = 1 and rk{\(j)')) > d. Then |0') has a unique 
representation 10') = j^b , for some complex 

numbers a,;, 1 < i < n, with w = J27=i I'^iP — ^- 
\P) ~ J27=i^i\^i)c/ If Charlie applies the projec- 
tive measurement {\P) {P\, Ic — \P) {P\} , with probabihty 
1/w > he observes \P) and Alice and Bob are left with 
10'). Alice and Bob can then convert |0') into |0) with 
nonzero probability because the target state's Schmidt 
rank is not higher ^ . □ 



We note that the above result generalizes the already 
established condition of SLOCC convertibility between 
bipartite pure states of rcf. . If the initial joint state 
of Alice and Bob is pure then the transformation becomes 

H')ab\0)c ^^'^ |0)ab- In this case, \i')AB is the only 
state in supp{pAB) and the transformation is possible if 
and only if rfc(|'0)AB) > rk[\(j))AB)- 

A unidirectional protocol like that described above is 
often called "one-shot" as Charlie's involvement consists 
of just making a measurement with rank-one measure 
operators and broadcasting the result to Alice and Bob. 
Hence \ip) abc can be converted to \4>) ab with a nonzero 
probability if and only if it can be done so by a one- 
shot protocol. The situation is strikingly different in the 
case of deterministic transformations since there exist tri- 
partite to bipartite conversions that require bidirectional 
collaboration between the parties in order to occur with 
probability one Q . 

According to Theorem 1, the problem of deciding con- 
version is reduced to whether a Schmidt rank d state 
exists in some subspace of Ha ® Hb- This question is a 
generalization of one sometimes referred to as Edmonds' 
Problem 0]: if M{dA, rfs) is the linear space of x dA 
matrices with complex coefficients and V is some sub- 
space of AI{dA, dB), decide whether there exists a rank 
d = Tahi{dA, dB} (i-e. full rank) matrix in V. Intuitively, 
one would expect that an upper bound exists on the di- 
mension s of subspaces containing only states of rank 
strictly less than d. Indeed, Flanders provides the bound 
s < d ■ maic{dA, ds} EEll- Thus, 



Corollary. // rfc(|0)) = d and dim[supp{p 
max{dA,rfB}, then \iP)abc ""^^^ 14') ab- 



4, X 

ab) 



> d 



It requires more work when dim[supp(p^^)] < d ■ 
maxjd^jrfs}. As recognized by previous investigators, 
determining whether a matrix subspace is singular can be 
cast into a polynomial identity testing question 12, 13| . 
We generalize their approach to the subject at hand. Let- 
ting {Tli}i=i---n denote the nonzero subnormalized eigen- 
states of p^g in matrix form, any state in supp{p^g) can 
be expressed as 



n(u) = uilii 



(2) 



where u = (ui, • • • , u„) is an rt-dimensional complex vec- 
tor. Construct the following real-valued function 



5(u) = ^|detM|^ 



(3) 



where Kd ranges over the set of d x d sub-matrices of n(u) 
and det(Kd) denotes the determinant of k^- Note that g 
is a nonnegative polynomial of degree no greater than 
2d in the real variables {a^, &i}i=i...n where Ui = Oi + ibi. 
Then deciding whether g is identically zero is the same as 
determining whether a Schmidt rank d state is obtainable 
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from \'iP)abc since ^(u) = if and only if no matrix of 
rank at least d exists in the span of {ni}i=i...„. 

As mentioned above, Polynomial Identity Testing is 
a classic problem in theoretical computer science with 
many important applications In general, given 

two polynomials /(x) and p(x), it can always be decided 
if / = P by multiplying out the polynomials and check- 
ing whether their coefficients match. However, the num- 
ber of multiplications required for this procedure scales 
exponentially in the degree of the polynomials and at 
the present no sub-exponential deterministic algorithm 
is known for polynomial identity testing [13, [3] ■ 

On the other hand, if one relaxes the determinis- 
tic condition, randomized polynomial-time algorithms 
exist _that can decide with a high probability of suc- 
[l6l |. A standard algorithm uses the Schwartz- 
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cess 

Zippcl lemma which states that for some n-variate poly- 
nomial f{xi, ■ ■ ■ ,Xn) over a field K and having degree 
no greater than d, if / is not identically zero, then 
PToh[f{x[,--- ,x'j^) = 0] < where each a;'- is inde- 
pendently sampled from some finite set X C K. Thus, to 
test with success probability at least 1 — whether a 
polynomial / is identically zero, one evaluates / on values 
chosen from set X and decides a zero identity if and only 
it f{x'i, - ■ ■ ,x[^) — 0. Using this algorithm on the analy- 
sis of g in Eq. Q allows for the following classification 
of the tripartite to bipartite conversion problem. 

Theorem 2. There exists a polynomial time randomized 
algorithm that, given states \'iI')abc md \<j))AB, decides 
correctly if ABC ^^^'^^ \<t>)AB is feasible with probability 
> 2/3. 

In other words, the problem of deciding tripartite to 
bipartite SLOCC convertibility belongs to the complex- 
ity class BPP, which consists of decision problems solv- 
able by Bounded-error Probabilistic Polynomial time al- 
gorithms. The error probability 1/3 can be made ex- 
ponentially small by repeating the algorithm and out- 
putting the majority of the outputs of each repetition. 

The randomized algorithm described above can be 
used to construct an LOCC protocol that completes 
any feasible tripartite to bipartite transformation with 
nonzero probability. Charlie makes 2n independent 
samplings ai, 6i, a2, 62, • • • ,an,bn from the integer set 
{1, • • • , M}, where M is an integer larger than 2d. Then 
he constructs u = (ai -I- i&i,- • • , a„ -I- i6„) and evalu- 
ates g{u) in Eq. ([3]). If ^(u) is nonzero, the state 
|P) = -^Y^^=i^Mi) is formed with TV being the ap- 
propriate normalization factor. By the Schwartz-Zippel 
lemma, such a u can be found with success probability at 
least 1 — |j in the case that g is not identically zero, and 
an appropriate M can be chosen to make this probabil- 
ity sufficiently large. Then Charlie performs a projective 
measurement {\P){P\,Ic — \P){P\}^ and Alice and Bob 
will share the unnormalized state 11^/ (u) = X]r=i "^i^i 



with nonzero probability when the outcome is \P). They 
will then be able to probabilistically obtain the desired 
target state as is with Schmidt rank at least d. 
One drawback of the above procedure is that we need to 
evaluate <?(u). However, the explicit form of g is unknown 
and may be very complicated as we need to sum over 
all determinants of c? x d sub- matrices of H(u). Fortu- 
nately, we can avoid evaluating ^(u) directly by checking 
whether the matrix rank of H(u) is larger than d, which 
can be done efficiently in polynomial time of n and dA( b) ■ 
It is interesting how the problem can be turned around 
by reducing any polynomial identity testing question to 
a decision of SLOCC convertibility. In an important 
work by Valiant [l3| (slightly improved in he shows 
that any polynomial p{xi, • • • , Xm) over a field K of for- 
mula size e can be expressed as the determinant of some 
(e -I- 2) X (e -I- 2) matrix Hp(a;i, • • • ,x,n)with. entries in 
{xi, Xm} and the underlying field [I4I, and the con- 
struction of Hp from p can be done in polynomial time. 
Let Up = Uo + T,'^lX^n„ where H^ e ]^{e+2)x{e+2) ^ 

< i < m. We claim that span{Hi : < i < m} con- 
tains a nonsingular matrix if and only if Hq -|- spanjHj : 

1 < i < m} does. The "only if" direction is triv- 
ial. Assume that A = X^IHo '^i^i is nonsingular. If 
ao 7^ 0, then Hq -I- Xil^i i^ nonsingular. Other- 
wise, for sufficiently large k, det(Ho -I- kA) > 0. Thus 
Hq + span{Hi : 1 < i < m} contains a nonsingular ma- 
trix. 

To conclude the reduction, for each i, < i < m, 
let \ei) be the bipartite state corresponding to H,, and 
Pab = ^ 1]"=! Letting IV') ABC be a purifi- 

cation of Pab £md a-ny rank n bipartite state, 

p{xi,--- ,Xm) is not identically zero if and only if 

\iP)aBC ^^"^ \(f>)AB- 

In conclusion, the results of this Letter help contribute 
to the complexity hierarchy of SLOCC pure state trans- 
formations. For bipartite transformations, the ques- 
tion of convertibility reduces to matrix rank calculations, 
which can be done in deterministic polynomial time. As 
seen here, for tripartite to bipartite conversions, deter- 
mining feasibility is equivalent to testing the identity of 
a given polynomial, which can be done in randomized 
polynomial time, and whether or not the algorithm can 
be derandomized is a major open problem in theoretical 
computer science. For general tripartite transformations, 
no polynomial time algorithm (deterministic or random- 
ized) is known, and under the common belief in complex- 
ity theory, no polynomial time algorithm exists. 

The next question might be how the complexity of 
calculating optimal conversion probabilities compares for 
transformations involving a different number of parties. 
For bipartite pure state conversions, the problem is al- 
ready known to have a polynomial-time solution It 
is interesting how in all these cases, abstract computa- 
tional questions can be used to solve a seemingly unre- 
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lated physical problem of entanglement conversion and 
vice versa. Such results demonstrate the value of com- 
plexity theory in the study of quantum systems and re- 
veals an important relationship between the two fields. 
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